Magnon dispersion in the field-induced magnetically ordered phase of TICUCI3 
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The magnetic properties of the interacting dimer system TICUCI3 are investigated within a bond- 
operator formulation. The observed field-induced staggered magnetic order perpendicular to the 
field is described as a Bose condensation of magnons which are linear combinations of dimer singlet 
and triplet modes. This technique accounts for the magnetization curve and for the field dependence 
of the magnon dispersion curves observed by high-field neutron scattering measurements. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.40.Cx 
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TICUCI3 is an insulating quantum spin system with 
a gap in the spin excitation spectrum [Q at zero field 
which originates from dimerization of the 5=1/2 spins 
of the Cu 2+ ions. The compound is isostructural to 
KCUCI3 and the crystal structure can be considered as 
coupled two-leg ladders separated by T1+ ions. How- 
ever, inelastic neutron scattering (INS) measurements of 
triplet magnon excitations found that the magnon modes 
have significant dispersion in all three spatial dimen- 
sions for both TlCuCl 3 and KCuCl 3 [§ § [§. This indi- 
cates that these compounds are three-dimensional (3d) 
interacting dimer systems with interladder interactions 
stronger than the interdimer interactions within the lad- 
ders. As observed in INS experiments by Cavadini et 
al, the magnon modes are split into three by a mag- 
netic field, with splitting proportional to the field, and 
the lowest mode becomes soft at a critical field H c [||. 
For TICUCI3 (KCUCI3), the zero-field excitation gap and 
H c are respectively OJmeV (2.6meV) and 5.7T (20T). 

For H > H c , both compounds show a uniform magne- 
tization parallel to the field H. Elastic neutron scatter- 
ing measurements on TICUCI3 show in addition a stag- 
gered magnetic order perpendicular to the field 0. A 
Goldstone mode is then expected due to the breaking of 
rotational symmetry around the field axis. Very recently, 
Ruegg et al. observed such a gapless mode for H > H c 
in TICUCI3 ||, and also reported a renormalized field 
dependence of the higher magnon modes [^) . 

Field-induced magnetic order in otherwise gapped lad- 
der systems has been described jlO], [ll| by theoretical ap- 
proaches which focus only on the singlet and the lowest 
triplet magnon. Giamarchi and Tsvelik jt2) cast their 
theory as a Bose condensation of a soft mode, and for 
TICUCI3 the Bose condensation of magnons has been 
used Jl3|, [li) to account for the observed temperature de- 
pendence of the magnetization. A quantum Monte Carlo 
simulation on a simplified 3d cubic lattice was also in 
agreement with an interpretation as magnon condensa- 
tion fnj. 

We note that the structure of the observed staggered 
order is related to the wavevector of the soft magnon 
[^|, , in support of the idea of a magnon Bose conden- 



sation. In this paper we develop a microscopic theory 
of field-induced magnetic order which takes into account 
the higher triplet modes. These modes, as we will show 
below, must be included to obtain a complete description 
of the condensate and of the evolution of the magnon dis- 
persion in the presence of field-induced magnetic order. 
We use a bond-operator formulation which retains all 
four states of each dimer, analogous to the treatment of 
bilayer antiferromagnets by Sommer et al. Jl6| . 

For the parameterization of the couplings between the 
Cu 2+ ion spins, we follow Ref. 13]. The unit cell contains 
two equivalent dimers, and the Hamiltonian is given by 

j 

+ Ji Sn,j ■ S n,j+d 1 + J 'l S lj ■ Srj+di 

Jsisl • + slj ■ sl j+d3 _)] + [1 <- 2]. 



Here S™^ is the spin 5 = 1/2 operator in unit cell j on 
the sublattice m = 1,2, and n( — l,r) denotes the left 
or right spin of the dimer. d\ — a, = 2a + c, and 
d$± = 6/2 ± (a + c/2), where a, b, and c are unit vectors 
corresponding to the a, b, and c axes, respectively ||. 

Because the intradimer exchange coupling J is the 
largest [|[ ||, we introduce bond operators and tj, 
based on these dimers [|l7|. In the presence of an ex- 
ternal field, the appropriate operators are |18| 

s f |o) - (IU)-Ut»/A 4|o) = -|TT), 
4|o) - (I U) + 1 IT»A/2, tl|o) = |u>- 

These have Bose statistics and are subject to the con- 
straint s*s + ^2 a t^ta — 1 at each dimer {j, m}, where 
a = +,0,—. We introduce the Fourier transforma- 



tion C ; i t = (l/A0£ fcC ™T e 
where 2N is the total dimer number. 

We restrict our considerations to zero temperature, 
and begin with the low-field region (H < H c ). Here the 
dimer singlets have the lowest energy, so the s-bosons 



for the dimer operators, 



2 



are taken to be condensed. The s-operators are replaced 
by a c-number = so5 ko an d the local constraint 

relaxed to the global one s* s = N- J2 k , a ^a^ka Re " 
taining quadratic terms in the triplet operators, 

m,fc a 

+ E (^fc^ m fca +rf Q ^ fca +H.c.) 

+ (/fc^-fco + Sfc*S)*-*o + H - c -)/2}, (2) 

where / fc = J a cosk x + J a2c cos(2k x + k z ), g k = 
2J a b c cos(k x + k z /2) cos(k y /2), a — —a, and m = 2, 1 
for rn = 1,2. The effective interdimer interactions || 
are given by J a = J x - J{/2, J q2c = -J2/2, and 
>/abc = (J3 — Ja)/2i where we note changes in the signs of 
the parameters J a 2c and, for J[ > 2Ji, J a . Introducing 
the operators = (t^, ± *fe„)/v / 2 gives two indepen- 
dent (±) modes |20| whose eigenvalues are obtained by 
the Bogoliubov transformation 



= v / (^ + /fc±3fc) 2 -(/fc±.9fc) 2 - a 9l x B H. (3) 

The Brillouin zone lies between — 7r and 7r (unit lattice 
spacing) in each direction, and the two branches corre- 
spond to the two-sublattice system. We treat only the + 
mode in the expanded Brillouin zone (— 27r < k z < 2tt) 
in the z-direction, because the magnon dispersions obey 



the relations E 



fc±(0,0,2ir),a' 



We extract the ef- 



fective interactions (J, J a , J a 2c, Jabc) from the data at 
zero field for both TlCuCl 3 and KCuCl 3 (Fig. |l|), and 
these are listed in Table I. The values arc consistent with 
the results of Ref . Q] . They are also similar to those of 
Ref . H , but not identical because the current treatment 
[leading to (|J)] goes beyond a pure dimer description. 
The signs are consistent with the expectation that all 
intcrsite parameters {J} in (1) are antiferromagnetic. 

The three magnon modes are degenerate in zero field, 
and, as shown in Fig. are split linearly by an external 
field in agreement with INS results §. Below H c , a t\_ 
triplet excited from the singlet condensate may propagate 
due to the interaction between triplets and singlets. The 
wave function of this excited state can be approximated 
by a linear combination of dimer singlets and triplets as 



«(| U)-UT))-«e 



i{k-rf-E k t) 



+ IJ I TT>, (4) 



where u is of order unity and v is a small, real coefficient. 
The expectation values of the spin operator components 
at a given dimer j, to in the state (jj) are 



(Sl.y)j. m 



~ {S r ,y)j,m 



- {uv/2)cos(k-rf 
~ —(uv/2) sin(fc • r r : 



(5) 

- E k+ t). 



This excited mode has a very small, uniform magnetic 
moment parallel to the field, and thus gains Zeeman en- 
ergy. Perpendicular to the field, it also possesses a finite 




FIG. 1: Zero-field magnon dispersion relations for 

TICUCI3 and KCUCI3. The rz;-axis labels represent the 
reciprocal-space points B=(0,27r,27r), C=(0,0,27r), D=(0,0,0), 
A=(7r,0,0), and F=(7r,0,27r). Experimental results [| are 
shown as points, and were measured at T = 1.5K (5K) for 
TlCuCb (KCUCI3). Solid lines are theoretical results using 
the parameters listed in Table I. 



magnetic moment, which is characterized by the wave 
vector k and energy E^ + , and is staggered (spins I and 
r oppositely aligned) due to the antiferromagnetic in- 
tradimcr coupling. For the mode t_, the direction of 
the uniform magnetic moment is antiparallel to the field, 
leading to a higher Zeeman energy. Finally, the t mode 
has no magnetization perpendicular to the field, and a 
moment parallel to the field which is modulated with 
wavevector fc, so its energy does not shift with the field. 
On increasing the field, these modes shift position with- 
out changing the shape of their dispersion (|j|), and the 
lowest (a = +) mode becomes gapless at the point C 
with Q = (0, 0, 27r) [Fig. 0(a)], which determines the crit- 
ical field. Thus at H = H c , the lowest mode exhibits a 
quadratic dependence on fc around Q. 

To describe the field regime with H > H c , the Hamil- 
tonian ^ must be extended to include triplet-triplet in- 
teractions. The coefficients of these terms involve combi- 
nations of the intersite exchange constants in (1) beyond 
the interdimer interactions in (0). To determine these we 



TABLE I: Effective dimer interactions. 





TlCuCl 3 


KCuCl 3 


J 


[meV] 


5.501 


4.221 


Ja 


[meV] 


-0.215 


-0.212 


Jalc 


[meV] 


-1.581 


-0.395 


Jabc 


[meV] 


0.455 


0.352 
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FIG. 2: Calculated magnon dispersions for TICUCI3 at H = 
H c (= 5.6T) (a), and H > H c (H = 12T) (b). 

have made the simplifying assumption that J\ = J3 = 0, 
so that the three remaining coefficients, J[, J' 2 , and J3, 
are specified by the interdimer interactions. Because 
the additional triplet-triplet interactions are largely gov- 
erned by terms involving J' 2 , this is not a significant ap- 
proximation. A further assumption is made by neglecting 
terms involving three t operators, which give only small 
corrections concentrated in the region of maximum stag- 
gered magnetic order. For H > H c , Bose condensation 
of the lowest triplet implies a macroscopic occupation of 
the ty, + mode at Q. 

For a full description of this regime, the singlet and 
triplet operators are transformed |16| to 

< = us k+^ +Q++ yt^ +Q _), 
b% + =u(xt% + +yt%_)-vs% +Q , (6) 

j m _ ,m h m ~ Tt m — iif m 

°kQ - r fcO' °fe- — xl k- V l k+- 

The fc-independent coefficients u, v, x, and y arise from 
two unitary transformations, and may be written as 
u = cos 9, v = sin#, x = cos (ft, y — sin0, with 9 and 
4> to be determined. We treat the operator as uni- 
formly condensed, = ao<5fccb and the ground state as 
a coherent condensate of the a™ operator. We empha- 
size that the highest triplet mode (£q_ ) also participates 

in the condensate, because Tt contains processes t ,t'_ss 
which nucleate i+ and t_ triplets from singlets. The lin- 
ear combination of singlet and triplets in the condensate 
a™ yields a staggered magnetization perpendicular to the 
field with wave vector Q, as observed in Ref. 0. 

With the transformation of Eq. (]?]), the Hamiltonian 
to quadratic order in the b operators takes the form 
H = 0(b°) + OO 1 ) + 0(b 2 ). Because the particle num- 
ber is unaltered, the c-number ao may be replaced by 
a* a = N — J2k, a b kJ b ka The constant terms (0(6°)) 
represent the mean-field energy of the ao-condensate, 
and the parameters (0, (j)) are chosen to minimize this en- 



FIG. 3: Field dependence of the energy gap of the three 
magnon modes in TICUCI3. Points are data from the INS 
experiment of Ref. || , and solid lines the results of the present 
theory. Inset: field dependence at high fields. We have used 
g = 2.16 for the g-factor of TlCuCFs |§. 

ergy, which also eliminates the 0(b l ) terms in the trans- 
formed Hamiltonian. The critical field H c is determined 
by the condition 9 — * 0, which gives purely singlet char- 
acter to the condensate. The limit (9 = ir/2, = 0) 
gives a condensate with purely triplet t + character, and 
determines the saturation field H s for full parallel polar- 
ization. The values for H c and H s are 

H c = VJ 2 ~ 2.7(1 J« I + \Ja2c\ + 2J abc )/(g» B ), 

H s = (J + 2|J a |+2|J a2c |+4J Qbc )/( 5MB ), (7) 

where the form of H c coincides with the soft-mode con- 
dition in the low-field regime. 

The 0(b 2 ) terms are diagonalized by a Bogoliubov 
transformation which yields the energies of the collec- 
tive modes of the condensate. Above H c , the lowest 
mode remains gapless, but develops a linear dependence 
on k near Q [Fig. ||(b)]. This is a Goldstone mode: a 
staggered magnetic moment M± , whose mean-field value 
is M± — uv{x + y)/\/2, is induced perpendicular to 
the magnetic field, and this breaks rotational symmetry 
around the axis parallel to the field. Rotations of this 
induced staggered moment are realized by changing the 
phase of x and y in Eq. © {x -> e-*x x,y —> e lx y, where 
X is the rotation angle), and do not change the energy. 
As a result the Goldstone mode remains gapless, and the 
field dependence of the higher modes is also renormalized, 
in agreement with experiment (Fig. ^|). The energy gaps 
of the higher modes Egp and E g ^ show an abrupt increase 
in slope at H c (Fig. [3]). At the saturation field H s , the 
k dependence of the lowest excitation mode near Q be- 
comes quadratic again. In the high-field region above H s , 
the condensate consists only of the lowest-lying triplet, 
and a gap reopens in the spectrum of the lowest (pure 
singlet) excitation mode (inset Fig. ||). 
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FIG. 4: Normalized magnetization curves for KCuCLj (a), 
and TlCuCi3 (b). The solid line (My) is the magnetization 
parallel to the field and the dashed line (Mj_) the square of the 
magnetization perpendicular to the field. Open (filled) points 
are experimental data for Mm (Mj_) measured at T — 1.3K 
(0.2K) @, ||, where g = 2.29 is the g-factor for KCuCl 3 . For 
TICUCI3 we have used the same g as in Fig. |§| 



We may also consider the parallel and perpendicular 
magnetization curves. For TICUCI3 (KCUCI3), the pa- 
rameters of Table I give a critical field of H c = 5.6T 
(19. 6T), which is consistent with the measured value 
For both TlCuCkj and KCuCl 3 , the square of the stag- 
gered moment (M?) has linear dependence close to H c 
and to H s (Fig. |j, indicating that M± ~ y/H - H c . 
For KCUCI3, the magnetization parallel to the held (Mm) 
is almost linear in H, as shown in Fig. |](a), in good 
agreement with experiment. For TICUCI3, the field de- 
pendence of M|| appears not to be linear in H, again 
in very good agreement with the observed form |J. The 
theoretical mean-field value is given by M\\ = v 2 (x 2 — y 2 ), 
from which it is clear that this difference is due to the 
magnitude of the interdimer interactions. If the contri- 
bution to the a™ operator of the highest triplet mode 
is neglected [Eq. ([?])], Mm would become completely 
linear in H . However, the strong interdimer interactions 
involve the highest mode largely through interactions of 
the type t+t^ss, which cost Zeeman energy, thus reduc- 
ing the value of v for the a™ operator, and consequently 
M|| is suppressed near H c for TICUCI3. 

We have studied the evolution of the magnon disper- 
sion in TICUCI3 and KCUCI3 as the magnetic field is 
tuned through the quantum critical point separating a 
gapped spin-liquid state from a state of field-induced 
staggered magnetic order which exists between the crit- 
ical and saturation fields. A consistent theory requires 
the admixture of both the lowest and highest triplet 
modes into the singlet dimer state to form a Bose con- 
densate. The spectrum has a gapless Goldstone mode 
associated with the breaking of rotational symmetry by 



the staggered magnetic order, as observed by Riiegg et 
al. H for TICUCI3 . The two higher excitation modes are 
also renormalized, in further agreement with observation 
Finally, our zero-temperature mean-field descrip- 
tion is also in good accord with measurements || |?]. Q 
of the uniform and staggered magnetization for both 
KCuChj and TlCuCl 3 . 
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